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ABSTRACT 


This  report  documents  the  results  of  a  test  program  on  a  Rotating  Pendulum  Accel¬ 
erometer  (RPA)  to  determine  its  applicability  to  the  Dormant  Missile  Inertial  Guidance 
and  Control  System. 

The  RPA  has  attracted  attention  in  this  program  because  its  accuracy  is  independ¬ 
ent  of  elastic  properties  of  matter,  critical  geometry,  and  slowly  varying  parasitic 
forces.  Accuracy  is  dependent  only  upon  the  repeatability  of  an  angular  null  Indication, 
and  the  ability  to  measure  time. 

Tests  were  performed  primarily  to  determine  stability  and  quick  starting  charac¬ 
teristics.  The  feasibility  model  originally  fabricated  was  modified  to  attain  the  degree  of 
quality  expected  of  future  models.  The  tests  demonstrated  feasibility,  and  performed  in 
the  manner  predicted,  showing  the  following  characteristics: 

-  Low  noise 

-  Stability  required  by  dormant  system 

-  Ability  to  start  quickly 

-  Capability  of  perfect  temperature  compensation 

On4he  basis  of  these  conclusions,  it  is  strongly  recommended  that  the  following 
program  be  authorized: 

-  Design  and  fabrication  of  an  engineering  prototype  RPA  (preferably  with  a 
synchronous  drive  to  simplify  *'->.mputation,  etc.) 

-  Functional  and  environmental  test  of  the  engineering  prototype 

-  Study  of  the  possibility  of  modifications  to  the  acceleration  sensor,  auxiliary 
equipment,  and  data  processing  methods.  It  is  intended  that  these  modifications 
would  lead  to  simplification  of  equipment,  improved  performance  and  accuracy, 
and  shorter  reaction  time. 
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SECTION  1 
INTRODUCTION 


This  report  presents  the  results  of  a  test  program  to  determine  the  applicability 
of  the  Rotating  Pendulum  Accelerometer  (RPA)  to  the  requirements  of  the  Dormant  Missile 
Intcrtlal  Guldanceand  Control  System.  The  RPA  Is  particularly  attractive  for  the  dormant 
missile  application  because  its  accuracy  is  independent  of  the  elastic  properties  of 
mntter,  critical  geometry,  or  slowly  varying  parasitic  forces.  The  novel  approach  of 
the  RPA  makes  Its  fundamental  accuracy  dependent  only  upon  the  repeatability  of  an 
angular  null  Indication  and  the  ability  to  measure  time. 

1.1  BACKGROUND 

The  RPA  was  conceived  by  Ford  Instrument  Co.  In  1958  and  was  patented  (No.  2, 
030,624)  in  May  1900,  Initial  Interest  was  retarded  to  some  extent  because  of  the 
novelty  of  the  npproach  and  its  dependency  upon  high  digital  counting  rates.  However, 
over  the  past  few  years,  the  state-of-the-art  of  digital  circuitry  has  advanced  to  the 
point  where  circuitry  capable  of  counting  at  the  required  rates  Is  readily  available. 

This  has  resulted  In  renewed  interest  in  the  RPA.  In  May  1902,  Ford  Instrument  Co. 
Initiated  activity  to  fabricate  and  test  a  feasibility  model  of  the  RPA.  >The  results  of 
the  activity  clearly  demonstrated  the  feasibility  of  the  RPA,  and  verified  the  analytical 
prediction  regarding  performance,  within  the  limits  of  accuracy  of  the  test  equipment. 

1.2  OBJECTIVE  AND  SCOPE 

Thc-objectlve  of  the  subject  RP,4  Jest  program  was  to  evaluate  the  application 
of  the  RPA  to  the  dormant  missile  requirements.  This  was  to  be  accomplished  by 
performing  tests  on  the  feasibility  model  to  Investigate  Its  stability  and  quick  start 
cluiracterlstlcs.  Inasmuch  as  the  RPA  feasibility  model  was  not  intended  to  demon¬ 
strate  the  accuracy  potential  of  the  RPA,  the  accuracy  of  the  feasibility  model  cannot 
be  considered  as  representative  of  what  can  l>e  obtained  In  a  prototype  unit.  However, 

In  order  to  investigate  adequately  the  stability  and  quick  start  characteristics  of  tin* 

RPA,  several  modifications  were  made  (at  no  cost  to  the  program)  to  both  the  feas¬ 
ibility  model  and  the  associated  test  equipment.  In  addition,  an  analysis  ol  several 
possible  sources  of  error  was  performed  (at  no  cost  to  the  program)  and  Is  reported 
in  Section  5. 


SECTION  2 

BASIC  CONSIDERATIONS 


The  basic  considerations  involved  in  the  application  of  the  RPA  are  discussed  in 
the  following  paragraphs. 

2.1  PRINCIPLE  OF  OPERATION 

The  Rotating  Pendulum  Accelerometer  represents  a  radical  departure  from  the 
conventional  approach  to  the  measurement  of  acceleration.  Whereas  all  known  existing 
accelerometers  involve  a  balance  of  forces  to  constrain  the  seismic  mass  from  moving, 
the  RPA  approach  allows  an  applied  acceleration  to  influence  the  motion  of  a  seismic 
mass,  and  computes  acceleration  from  measurements  of  the  motion.  The  seismic  mass 
used  in  the  RPA  is  an  unbalanced  cylinder  rotating  on  a  hydrostatic  gas  bearing  which 
represents  a  close  approximation  to  an  ideal  frictionless  bearing.  The  effect  of  an 
acceleration-normal  to  the  axis  of  rotation  is  to  modulate  the  angular  velocity  of  the 
rotating  cylinder.  This  modulation  is  detected  in  the  RPA  by  measuring  the  time  re¬ 
quired  to  rotate  through  successive  quarter  cycles  of  rotation. 

2.2  PERFORMANCE  EQUATION 


It  may  be  shown  that  if  the  applied  acceleration  is  constant  over  a  cycle  of 
rotation,  then  a  measurement  of  the  time  required  to  traverse  quarter  cycles  is 
sufficient  to  determine  the  components  of  acceleration  along  two  reference  axes  normal 
to  the  axis  of  rotation.  To  a  first  approximation,  the  RPA  behaves  like  a  linear  device. 
Thus,  the  acceleration  components  mi •>  be  computed  from  the  following  equation. 


where  A  ■  Acceleration  component 
x  3 

K  ■  Constant  =  2  r  PE 


K 


At 


o>  =  The  natural  radian  frequency  of  the  RPA  when  oscillating  through  smnll 

r 

angles  in  the  reference  gravitation  field  gg 

At  »  Difference  in  time  required  to  traverse  successive  half  cycles  of 
rotation. 

T  ■  Period  of  rotation. 
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The  error  involved  in  the  use  of  this  equation  depends  upon  the  magnitude  of  the 
modulation.  Using  At /r  as  a  modulation  index,  the  equation  has  an  error  of  approxi- 
mately  1  part  in  10  when  At/ T  is  less  than  0.03.  More  exact  equations  yielding  any 
desired  accuracy  may  be  derived  by  adding  higher  order  terms  to  the  linear  approxi¬ 
mation.  Based  upon  an  analytic  investigation  of  the  undamped  case,  an  expression  for 
acceleration  having  errors  of  about  3  parts  in  10  has  the  following  form: 


Although  this  equation  was  derived  from  consideration  of  the  undamped  motion  of 
the  RPA,  it  is  expected  that  a  solution  of  the  same  form  is  equally  applicable  to  the 
damped  case,  since  the  magnitude  of  the  damping  in  a  hydrostatic  gas  bearing  (as  shown 
in  paragraph  5.2)  has  a  very  small  effect  upon  the  magnitude  and  phase  of  the  modulation. 

2.3  ACCURACY 

The  results  of  an  error  analysis  of  the  RPA  indicate  that  the  main  sources  of 
error  are  parasitic  torques  which  vary  at  a  rate  comparable  to  the  rate  of  rotation, 
quantization  errors  due  to  measurement  of  finite  increments  in  time,  and  the  stability 
of  the  natural  frequency  of  the  pendulum  and  of  the  clock  used  to  measure  time.  Slow 
variations  in  the  torque  used  to  sustain  rotation,  changes  in  viscosity  of  the  gas  due  to 
temperature  changes,  and  changes  in  dampii  >:  due  to  changes  in  bearing  eccentricity 
under  load  (see  paragraph  5,2)  cause  negligible  error  in  the  computed  acceleration. 

In  a  temperature-compensated  RPA,  it  is  expected  that  changes  in  temperature  will 
also  cause  negligible  error  (see  paragraph  5.1). 

Since  the  RPA  yields  a  computation  of  acceleration  or  a  velocity  increment  for 
each  cycle  of  rotation,  the  rapidly  varying  parasitic  torques  and  quantization  errors 
cause  what  may  be  considered  noise,  i.c.,  random  variations  in  the  velocity  increment 
about  its  true  value.  Instability  in  the  pendulum  natural  frequency  or  clock  causes 
changes  in  scale  factor.  The  major  source  of  error  in  a  single  velocity  increment  is 
noise.  However,  when  a  large  velocity  Increment  is  measured  by  adding  many  small 
increments,  the  percentage  error  due  to  noise  varies  inversely  with  the  square  root  of 
the  number  of  increments,  due  to  the  partial  cancellation  of  error  that  occurs  when 
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increments  having  positive  errors  are  added  to  increments  having  negative  errors. 

After  a  very  large  number  of  increments,  the  error  is  determined  only  by  the  stability 
of  the  pendulum  natural  frequency  and  clock.  It  may  be  shown  that,  for  equal  increments, 
the  total  error  is  given  by: 


(eV)rms 

.jr 

(cAV)rms 

2 

V 

Vn 

AY 

where 


(eV|rms 


Normalized  rms  error  in  velocity 


n  =  Number  of  velocity  increments 

(eAV)rms  =  Normalized  rms  error  in  one  velocity  increment  due  to  noise 

aV 

K  =  Scale  factor  error. 


A  plot  of  this  equation  for  a  ten-to-one  ratio  of  noise  error  to  scale  factor  error 
per  increment,  for  equal  Increments,  is  shown  in  figure  2-1.  Based  upon  the  test 
results  to  date  and  state-of-the-art  of  digital  circuitry,  it  is  estimated  that  the  normal¬ 
ized  noiffS'would  be  about  5  x  10"5,  assuming  the  RPA  rotates  at  about  25  cps.  This 
figure  is  sufficiently  low  to  make  stability  the  determining  factor  with  regard  to  accuracy 
after  a  large  number  of  velocity  increments. 
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Figure  2«1.  NurmallrH  Velocity  Error  vermin 
NumNT  of  Increment* 


SECTION  3 

DESCRIPTION  OF  TEST  SET-UP 

3.1  PARAMETERS  OF  RPA 

The  parameters  of  the  RPA  used  In  this  test  program  were  as  follows: 

-  Natural  frequency:  1  cps 

-  Pendulosity:  67  gm-cm 

2 

-  Inertia:  1725  gm-cm 

-  Air  bearing  time  constant:  163  sec. 

3.2  DESCRIPTION  OF  TEST  EQUIPMENT 

The  test  equipment  used  in  the  performance  of  those  tests  was  considered 
adequate  to  establish  feasibility  of  the  RPA  for  the  dormant  application.  Figures  3-1 
and  3-2  are  a  sketch  and  a  photograph,  respectively,  of  the  RPA.  Figure  3-3  shows 
a  photograph  of  the  test  set-up.  A  brief  description  of  the  test  equipment  follows. 

A  stable  base  was  selected  on  which  a  sine  bar  was  placed.  The  accelerometer 
was  then  securely  fastened  to  the  sine  bar,  which  permitted  the  RPA  to  be  accurately 
positioned  with  respect  to  the  vertical,  through  the  use  of  calibrated  gage  blocks.  The 
half  cycle  outputs  of  the  accelerometer  were  then  inserted  into  two  Hewlett-Packard 

a 

time  interval  counters  (accurate  to  l  part  in  10  ),  which  measured  the  time  required 
to  travorsiTPShSecutlve  half  cycles.  These  data  were  then  recorded  by  digital  re¬ 
corders,  and  later  processed  nnd  analysed. 
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Figure  3-1.  Rotating  Pendulum  Accelerometer 
(Feasibility  Model),  Line  Drawing 


Figure  3-3,  nutating  Pendulum  Arrelrromrter 
(Feasibility  Model) 


SECTION  4 
TEST  RESULTS 


The  laboratory  test  program  conducted  using  the  previously  described  model  of  the 
RPA  and  its  associated  test  set-up  verified  the  excellent  stability  and  quick  start  charac- 
teristics-of  theflPA  within  the  limitations  of  the  test  set-up.  Due  to  certain  inadequacies 
in  the  physical  test  set-up  and  resolution  errors  in  the  clock  reference,  it  was  impossible 
to  demonstrate  performance  with  the  desired  accuracy.  Physical  variations  (to  be  de¬ 
scribed  in  paragraph  4.3)  caused  fluctuations  in  input  acceleration  that  were  not  antici¬ 
pated,  while  the  clock  resolution  error  caused  an  error  in  At,  varying  between  plus  and 
minus  2  microseconds.  Since  the  full-scale  value  of  At  (for  0.03  modulation  and  r  -  0.4 
sec)  was  12,000  microseconds,  the  error  limits  in  At  (and,  effectively,  in  At/r**)  were 
±1.7  x  10*4  due  to  clock  resolution  alone.  Averaging  of  successive  measurements  of 

At/r^  for  constant  input  acceleration  should  reduce  this  error  effect.  A  statistical 
data  analysis  was  performed  on  measured  data  at  two  RPA  speeds,  and  is  described  in 
paragraph  4.1 . 

Additional  tests,  described  in  this  section,  include  calibration  runs,  linearity  tests, 
quick  start  test,  and  stability  test.  The  calibration  runs  were  performed  in  order  to 
determine  the  table  bias  error,  the  bias  error  due  to  RPA  half  cycle  unbalance,  and  the 
linear  coefficient  of  the  RPA.  The  linearity  run  was  to  be  used  to  determine  the  coeffi¬ 
cients  of  higher  order  terms  in  (  At/r*)  to  correct  the  errors  in  the  linear  approxima¬ 
tion.  Unfortunately,  lack  of  time  made  it  impossible  to  include  these  data  in  this  report. 
However,  these  data  will  be  forwarded  .-'e  soon  as  they  become  available.  The  quick  start 
test  demonstrated  how  rapidly  the  RPA  output  settles  out  to  its  final  value  after  starting 
from  a  dormant  mode.  The  stability  test  demonstrated  the  repeatability  of  the  RPA  for 
a  constant  input  over  an  extended  period  of  time. 

4.1  PROBABILITY  DISTRIBUTION  OF  NOISE 

As  discussed  in  paragraph  2.3,  the  RPA  yields  a  computation  of  acceleration  or 
a  velocity  increment  for  each  cycle  of  rotation.  Each  of  the  velocity  increments  will  have 
random  noise  superimposed  upon  the  true  value  generated  by  rapidly  varying  parasitic 
torques,  and  by  the  quantisation  of  time  due  to  the  measurement  of  time  in  finite  Incre¬ 
ments  (one  microsecond  In  the  test  set-up).  \ 
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The  magnitude  of  the  noise  was  measured  at  rotation  rates  of  1  cps  and  2.5  cps. 
Based  on  100  consecutive  measurements,  the  probability  frequency  function  was  plotted 
with  normalized  (with  respect  to  the  1(T  value)  deviation  from  the  mean  as  abscissa.  The 
results  are  shown  in  figures  4-1  and  4-2.  As  may  be  seen  by  comparing  the  measured 
distribution  with  theoretical  values  for  the  normal  distribution,  the  noise  is  very  nearly 
normally  distributed.  The  normalized  standard  deviation  from  the  mean  was  0.74  x  10"4 
at  1  cps  and  0.72  x  10"4  cps  at  2.5  cps. 

The  maximum  error  in  A  t  due  to  quantization  of  time  is  2a,  where  a  is  the  recipro¬ 
cal  of  the  clock  frequency.  Its  standard  deviation  is  2  a  A/T ;  hence,  the  normalized 
standard  deviation  due  to  quantization  of  time  is: 

2a 

0.03  /3~  T 

where 

°NC  “  Normallzed  standard  deviation  due  to  quantization 

a  *  Reciprocal  at  clock  frequency 
T  ■  Pendulum  period  of  rotation 
Substitution  of  numbers  into  the  equation  yields 

<rNC  «  1.5  x  10-4  at  2.5  cps 

ffNC  "  nt  *  cpB 

It  may  be  seen  that  the  theoretical  noise  due  to  quantization  is  of  the  same  order  of 
magnitude  as  the  measured  noise.  Theoretically  the  noise  due  to  random  torques  In  the 
pendulum  should  decrease  at  the  higher  speeds  (as  shown  In  paragraph  5.3)  while  the 
noise  due  to  quantization  should  Increase.  At  some  speed  they  will  be  approximately 
equal,  Based  upon  the  measurements  made,  It  appears  that  this  speed  wns  in  the  rnnge 
of  1  to  2.5  cps-forthtf  feasibility  model  and  the  associated  test  set-up.  The  nmgnltude  of 
the  noise  could  be  reduced  by  Increasing  both  the  speed  of  rotation  and  the  clock  fre¬ 
quency,  However,  due  to  the  cancellation  of  noise  errors  when  a  large  number  of  velocity 
increments  are  summed,  the  noise  level  In  the  feasibility  model  is  almost  adequate  for 
the  dormant  missile  application.  Furthermore,  It  Is  reasonable  to  expect  that  it  could 
be  substantially  reduced  as  a  result  of  a  development  program. 
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4.2  CALIBRATION  RUNS 

The  purpose  of  the  calibration  runs  was  to  determine: 

-  -  The  scale  factor  (of  the  RPA  linear  term). 

-  The  table  bias  (error  in  knowledge  of  the  vertical). 

-  The  RPA  half  cycle  bias  (measurement  error  due  to  the  difference  between  180- 
degree  half  angles  of  the  RPA  due  to  pick  off  error). 

The  data  wore  taken  for  very  low  modulation  (  At/T<0.003),  in  order  to  stay  within  the 
highly  linear  region  of  the  RPA. 

In  this  same  region  of  very  small  angular  deviations  from  the  vertical,  sin  0  may 
be  replaced  by  0,  and  the  normalized  acceleration  a/g  is  exactly  equal  to  0,  as  measured 
from  the  true  vertical.  Also,  it  may  be  shown  that  the  half  cycle  bias  causes  an  error  in 
the  measured  value  of  At/r3  equal  to  Kb/T2,  where  is  a  constant  equal  to  8/tt  ,  8 
being  the  radian  angle  difference  from  180  degrees  for  each  half  cycle  rotation.  Thus, 
the  measured  value  of  At/r2  may  be  expressed  as 

Mr  =  m0  +  b 

T3 

where 

K. 

b  =  7T  '  m0o 

and  0Q  is  the  table  offset  bias  (error  in  knowledge  of  vertical).  The  combined  bias  term 
b  will  be  constant  If  T  is  constant  for  a'\  data  points.  If  values  of  m  and  b  are  determined 
for  two  independent  calibration  runs  at  different  T  values,  a  pair  of  simultaneous  equa¬ 
tions  may  be  solved  to  find  values  for  K.  and  9  . 

b  o 

Calibration  runs  were  taken  for  r  =  0.4  second  (2.5  rps)  and  T  =  1.0  second  (1 
rps).  These  curves  were  fitted  by  a  least  square  error  linear  curve-fitting  technique  to 
determine  m  and  b  for  each.  The  point-by-point  deviation  of  each  data  value  from  the 
llnenr  fit  was  normalized  with  respect  to  full  scale  (output  for  0.03  modulation).  The 
calibration  curve  nnd  the  corresponding  normalized  error  curve  at  2.5  rps  arc  presented 
in  figure  4-3.  Each  data  point  actually  consists  of  the  average  of  10  consecutive  meas¬ 
urements  with  the  same  9  value.  At  2.5  rps,  the  rms  deviation  from  the  curve  fit  was 
CT  =  1,15  x  10*^,  These  deviations  nre  due  mainly  to  clock  resolution  errors 
and  also  other  test  uncertainties  which  will  be  discussed  in  paragraph  4.3. 
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The  zero  angle  reference  for  each  of  the  calibration  runs  was  taken  as  the  point 
of  minimum  measured  At/r  3.  These  reference  angles  were  different  for  the  two  runs 
due  to  the  Kb/T  3  bias;  the  difference  between  the  reference  angles  had  to  be  used  to 
modify  the  equations  used  to  isolate  Kb  and  9o.  The  value  determined  for  Kb  (6.100  x 
10’4)  corresponds  to  a  half  angle  error  of  approximately  6  minutes  of  arc. 

4.3  LINEARITY  RUNS 

It  had  been  planned  to  take  linearity  data  over  the  range  of  modulation 
-0.03  £  At/r  HZ  +0,03  and  to  curve-fit  these  data  using  a  linear  term  In  At/  t3  and 
some  higher  order  correction  terms  to  give  a/g  to  an  accuracy  of  1  part  in  103.  (The 
linear  curve-fit  Is  theoretically  good  to  1  part  In  103  over  this  range.)  Unfortunately, 
certain  unanticipated  uncertainties  in  the  laboratory  set-up  caused  dispersions  in  the 
measured  data  which  rendered  this  curve-fitting  impossible.  The  main  sources  of  these 
uncertainties  are  the  following: 

-  Shifts  in  the  mounting  table  which  caused  small  changes  in  the  effective  null 
positions. 

\ 

-  Play  in  the  sine  bar  mounting  frame. 

-  Play  in  the  mounting  of  the  RPA  (detected  near  the  end  of  the  program) 

-  Clock  resolution  errors . 

Additional  sources  of  small  variations  in  RPA  output  include  minute  fluctuations  in 
torquer  voltage  and  air  bearing  turbine  toroues. 

In  spite  of  the  above  error  sources,  the  RPA  did  exhibit  linearity,  within  the 
modulation  limits  specified,  of  1  part  in  103.  Howevor,  the  deviation  from  linearity 
did  not  Indicate  any  trend  that  could  be  curve-fitted  with  the  higher  order  terms.  For 
example,  the  error  plot  indicated  apparent  shifts  in  the  vertical  during  a  linearity  run. 

Most  of  the  major  error  sources  can  bo  corrected  readily,  and  reasonably  good 
linearity  data  should  be  available  in  a  short  time. 

4.4  QUICK  START  TEST 

The  quick  start  test  demonstrates  the  manner  in  which  the  RPA,  from  n  cold  start, 
rapidly  clostwlffah  Us  final  value.  The  test  was  run  with  n  constant  acceleration  Input. 
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The  RPA  was  energized,  spun  up  to  its  nominal  rotation  rate,  and  allowed  to  stabilize. 

Data  worcxccordedat  intervals  during  the  first  35  minutes  following  the  energizing. 

The  results  are  given  in  figure  4-4,  where  the  deviation  from  the  35-minute  value 
of  output  (normalized  with  respect  to  that  value  of  output)  is  plotted  as  a  function  of  time. 
Within  4  minutes,  the  RPA  output  was  within  a  range  of  deviation  of  3  parts  in  104  from 
its  final  value. 

4.5  STABILITY  TEST 

The  purpose  of  the  stability  test  was  to  demonstrate  the  repeatability  of  the  RPA 
output  over  an  extended  period.  The  RPA  was  rigidly  mounted  with  a  fixed  input  accelera¬ 
tion,  and  periodically  energized.  After  approximately  30  minutes  operation,  its  output 
was  recorded. 

Thc_aornialized  deviation  of  the  RPA  output  (plotted  in  figure  4-5)  was  found  to 
remnin  within  about  1.5  parts  in  104  from  its  mean  value. 
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SECTION  5 

ANALYTICAL  INVESTIGATIONS 


An  analysis  of  three  problem  areas  relating  to  the  RPA  is  presented  below. 
5.1  TEMPERATURE  STABILITY 


The  scale  factor  of  the  RPA  varies  inversely  with  the  square  of  the  pendulum 
natural  frequency.  The  problem  of  temperature-compensating  the  RPA  is  therefore 
very  similar  to  that  of  temperature- compensating  a  clock.  In  the  following  paragraphs 
it  is  shown  that  perfect  compensation  is  feasible  in  a  practical  design. 


The  linear  expression  for  the  acceleration  sensed  by  the  RPA  is 

I  o)a3  At 
a  4  mr 


mr 


2  7T 


3  At 


The  scale  factor  of  the  device  is  then  proportioned  to 


where  I  is  equal  to  the  inertia  of  the  floated  solid  cylinder  plus  the  inertia  of  the 
pendulous  mass,  and  mr  is  the  pcndulosity  of  the  device.  Any  variation  in  the  value  of 
•jjjp-  will  naturally  cause  a  corresponding  variation  in  the  computed  acceleration.  The 
largest  variation  in  1~r--  will  be  caused  by  temperature  variations;  this  error  can  be 
eliminated  by  selecting  the  system  parameters  so  that  they  conform  to  the  following 
expression: 


where  V  =  Coefficient  of  thermal  expansion  for  the  floated  cylinder 

q  =  Coefficient  of  thermal  expansion  for  the  support  member 
of  the  pendulous  mass 

m  »  Pendulous  mass 
r  =  Torque  arm  of  the  pendulous  mass 
M  =  Mass  of  the  floated  cylinder 
R  =  Radius  of  the  floated  cylinder 
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Examination  of  this  equation  indicates  that  the  value  of  must  be  less  than 
> 

-H-  since  — «-  is  a  positive  quantity.  A  simple  calculation  was  performed  in  order  to 
MIT 

determine  the  practicality  of  the  above  solution. 

in  1  ■  0 

Assuming  a  realistic  value  of  ■  j-,  and  values  of  v  =  11  x  10  /°C  and  q 

24  x  lO’V’C,  which  are  respectively  the  coefficients  of  the  thermal  expansion  for 
steel  and  aluminum,  and  R  =1.5  inches,  a  value  of  r  =  0.9682  inch  is  obtained.  These 
results  are  obviously  in  the  area  of  practical  design,  indicating  that  the  RPA  can  be 
designed  to  operate  with  no  variation  in  scale  factor  over  the  temperature  range. 
Figure  5-1  illustrates  the  design  configuration  which  could  be  used  in  temperature- 
compensated  RPA  units.  The  pendulous  mass  would  consist  of  a  material  with  a 
higher  coefficient  of  thermal  expansion  than  the  floated  cylinder. 

A  computation  was  performed  on  the  existing  RPA  in  order  to  determine  the 
variation  in  scale  factor  over  the  temperature  range.  A  variation  of  0.001%/°  C  was 
computed.  This  variation  is,  of  course,  not  indicative  of  the  scale  factor  sensitivity 
that  may  be  expected  of  a  temperature-compensated  RPA,  since  the  present  unit  was 
fabricated  from  parts  already  in  existence. 

5.2  EFFECT  OF  CHANGES  IN  DAMPING  DUE  TO  HIGH  g  ENVIRONMENT 

The  effect  of  damping  in  the  RPA  gas  bearing  is  to  reduce  the  amplitude  of 
modulation  by  a  factor  x  and  shift  the  phase  of  the  modulation  by  an  angle  £  *  where 
x  and  £  are  given  by: 


where 

Assuming 


f  =  tan’ 


1 _ 

T  )2 
2  TTT  p' 

^  2ttt  p 


) 


T  =  Period  of  rotation 
tp  =  Pendulum  time  constant 

T  =  163  sec 
T  =  0.05  sec 

r 


♦"Rotating  Pendulum  Accelerometer",  ARS  Journal,  April,  1901. 


then  substitution  of  T  and  Tp  into  the  equations  for  x  and  *  yields 

x  =  1.69  x  lO-9 
£  =  9.9  arc-seconds 

From  the  above  calculations,  it  is  evident  that  the  effect  of  damping  on  the 
magnitude  of  acceleration  is  so  small  that  even  an  order  of  magnitude  increase  in 
damping  would  still  produce  a  negligible  effect.  However,  a  substantial  change  in 
damping  would  appreciably  affect  the  phase  shift  and  introduce  cross  coupling  by 
the  amount  of  the  change  in  the  phase  shift.  The  results  of  an  analytical  investigation 
into  this  source  of  error  due  to  bearing  eccentricity  caused  by  a  high  g  environment 
are  presented  below. 

It  may  be  shown  that  the  equation  for  the  damping  torque  in  an  eccentric  bearing 
is  given  by: 


where 


1.  /I  -  8/h 
Vi'+T/h 

(8/h)2 


1  -  (8/h) 
1  -  (8/h)2 


Damping  torque 
Viscosity  of  gas 
Angular  velocity  of  bearing 
Radius  of  bearing 
Nominal  gap 

Displacement  of  bearing  from  its  neutral  position. 


The  first  two  factors  in  the  equation  lor  T^  result  from  damping  around  the  circum 
ference  of  the  bearing  while  the  third  term  arises  from  damping  on  the  bearing  end 
plates. 


Assuming  ^  =  -j-  (a  typical  value)  the  normalized  damping  torque  T^  is  given 


1  I  8  tan 


-1  /I  -  8  h 


DN  "  2.25 7T 


tan  VTTTh  + 

1  -  (  8  /h)  2 


.  -i  ITTti 

tan  ,  / — - —-X-  ,  7 r 

Vl  -  (8/h)2  +  T~ 

1  -  (  8  /h)  2 


where  Tj^  equals  the  damping  torque  normalized  with  respect  to  the  damping  torque 
at  zero  eccentricity. 
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A  plot  of  the  percent  increase  in  damping  as  a  function  of  eccentricity  is  shown 
in  figure  5-2.  From  the  curve  it  may  be  seen  that,  for  an  eccentricity  of  0.5,  the  change 
in  phase  shiftjLA.f)J.s:~ 

A£  =  (0.13)  9.9  arc-seconds 
=  1.3  arc-seconds 

The  above  computation  indicates  the  cross  coupling  error  due  to  bearing  eccen¬ 
tricity  under  load  is  acceptably  small  for  bearing  eccentricities  less  than  0.5. 

5.3  EFFECT  OF  RANDOM  VARIATIONS  IN  THE  SUSTAINING  TORQUE 

If  the  torque  used  to  sustain  rotation  of  the  RPA  changes  from  its  normal  value 
and  assumes  a  new  constant  value,  the  period  of  rotation  will  change.  But  since  the 
period  of  rotation  is  measured  and  used  in  the  computation  of  acceleration,  there  is 
no  error.  However,  if  the  sustaining  torque  changes  during  a  period  of  rotation,  then 
there  will  be  an  error  in  the  computed  acceleration.  The  results  of  an  analytic  inves¬ 
tigation  into  the  errors  resulting  from  a  random  variation  in  the  sustaining  torque  is 
presented  below.  The  analysis  is  based  upon  an  investigation  of  the  effect  of  the 
random  torque  on  At,  assuming  no  acceleration  input  and  assuming  the  RPA  is  in  a 
linear  range  of  operation  and  the  effect  of  damping  may  be  neglected.  It  is  felt  that 
these  assumptions  do  not  introduce  any  significant  errors  into  the  result. 

In  view  of  the  preceding  assumptions,  the  random  variations  in  the  angle  of 
rotation  may  be  computed  from  the  equation: 


where 


I  G  =  (t) 

I  =  Moment  of  inertia 
G  =  Angular  acceleration  of  rotor 
M(t)  =  Random  torque 


Since  it  has  been  assumed  the  RPA  is  in  a  linear  range  of  operation  (the  mod¬ 
ulation  is  small),  the  principle  of  superposition  applies  and  the  random  variations  in 
O  may  be  superimposed  upon  the  constant  angular  velocity.  The  equation  for,  A  t  is 
approximately 


At  = 
At  = 


(r  + 


('t  \  /  r  ,  °r/2  ' 

u -)  *  (-£- +  -rr_) 


°T  2°T/2 

“  11 


°r/2 

n 
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where 


At  =  Difference  in  time  between  successive  half  cycles  of  rotation 
T  =  Period  of  rotation 
n  =  Nominal  angular  velocity  of  RPA 
0T  =  Random  variation  in  8  at  the  end  of  one  period 
®T/2  =  ^anclom  variation  in  6  at  the  end  of  one- half  period 
The  above  approximation  solves  one  of  the  major  difficulties  in  analytic  investigations 
of  RPA  performance,  specifically,  the  problem  of  inverting  the  solution,  i.e.,  solving 
for  time  as  a  function  of  8.  The  mean  squared  value  of  At  Is: 

eM  .  *[£  (  0r  -  28r/2>  2] 

eM  -  -ij  E  [«t2-  46t  ®  T/2  *  «T/2J 

where  E  ( A t^J  =  expected  value  of  At^ 

Assuming  zero  initial  conditions,  and  recalling  that  10  =  M(t),  the  equation  for  the 
random  variation  in  8  is  as  follows: 

1  A 

9  'T  Jo  dT/o  M(l)dt 

In  the  proccnn  ot  cvnluntlnn  E|At2].  tt  In  nPccHsnry  to  know  the  nutocnrretntion 
function  of  the  random  torque.  To  obtain  a  solution,  an  exponential  autocorrelation 
function,  which  Is  good  approximation  to  the  autocorrelation  function  in  a  wide  variety 
of  random  processes,  was  assumed.  Thus: 

E  | M(tj)  M(t2)|  =  (T2  e  'A|ll  '  4| 

where  E  j  M(tj)  Mftg)  j  =  Expected  value  of  the  product  of  the  torque  at  times 

tj  and  t2 

2 

< 7  =  Mean  squared  value  of  6 

X  «  Constant,  which  determines  how  fast  the  random  torque  is  changing 
Evaluation  of  E  |  At“*|  in  the  manner  previously  described  is  a  long  and  tedious  lank. 

The  result  Is: 

«.  I  <  \>\2|  (T  2  .  6  1  2  6  2e*AT  8e'Ar/2 

E  I  (At)  |  =  - —K — K~  T  7m  rT  *  - ?  - a  +  - T  *  - T 

(2  TT  )2  I  2  6(AT)  (XT)3  (XT)4  (XT)4  (A TT 

•J 
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The  normalized  acceleration  error  is: 


g  A 
m 


( At)rms 

r3 
"  6.0  3 
“r7 


(At)rms 

0.03r 


where  e  A  =  Acceleration  error  normalized  with  respect  to  maximum  acceleration 
An  - - input  nt  A  t/T  =  0.03 


( At)rms 


1 


A  t 


F  A 

Substituting  ( At) r ms  into  the  equation  lor  yields: 

m 


gA 


m 


a  t 

0.03  (2 tF)\  ~r 


^  ■  <7^7  +  (XT, 


96  32e‘AT  128e‘AT/2 


i  1/2 


(XT) 


(XT)1 


It  may  be  noted  that  for  a  given  value  of  XT  ,  the  above  equation  predicts  that  the 
normalized  error  will  vary  directly  with  the  square  of  the  period  of  rotation.  Hence, 
one  wouid  expect  better  performance  at  higher  angular  velocities,  everything  else 
being  equal.  This  trend  was  observed  in  the  one-  <T  values  of  the  measured  noise  at 
1  cps  and  2.6  cps,  However,  the  improvement  in  noise  caused  by  parasitic  torques 
at  2.5  cps  was  probably  masked  by  the  quantization  noise  which  could  nccount  for  the 
total  magnitude  of  the  noise  at  2.5  cps  (approximately  10"4). 

A  plot  of  the  relative  noise  (square  root  of  the  quantity  in  brackets)  as  a 
function  of  (XT  )  for  constant  T  is  shown  h  .figure  5-3.  By  expanding  the  exponential 
terms  in  a  power  series,  it  may  be  shown  that  a  good  approximation  for  XT  «  1  is 
given  by: 


W?2'7r)T 


1  *  "sx  +  ~w  (Xt)2 


It  may  be  noted  that  this  equation  predicts  that  the  normalized  error  will  reach 
a  constant  value  as  (X  T  )  approaches  zero,  corresponding  to  a  constant  torque.  This 
result  is  Incorrect  and  is  caused  by  the  fact  that  A  t  was  computed  during  the  first 
cycle  of  rotation  after  the  torque  was  applied.  The  correct  steady-state  result  should 
be  computed  by  performing  the  required  integrations  over  a  cycle  starting  from  t^ 
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and  then  allowing  t  to  approach  infinity.  It  is  expected  that  the  steady- state  solution 
would  show  a  peak  at  some  value  of  (  XT  )  and  approach  zero  as  (XT)  approaches  zero 
and  infinity.  It  seems  reasonable  to  expect  the  peak  to  occur  at  (  X  T  )  =  1,  i.e.,  when 
the  noise  is  changing  at  a  rate  comparable  to  the  rate  of  rotation.  The  anticipated 
peaking  of  the  relative  error  is  shown  dotted  in  figure  5-3. 
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SECTION  6 

CONCLUSIONS  AND  RECOMMENDATIONS 


6.1  CONCLUSIONS 

The  principal  results  of  the  RPA  test  program  are  listed  below. 

-  The  normalized  (with  respect  to  maximum  acceleration  input)  root  mean  square 
value  of  noise  in  a  single  velocity  Increment  was  0.7  x  10‘4.  The  noise  was 
normally  distributed  about  the  mean  and  had  a  magnitude  very  nearly  equal  to 
that  which  results  from  the  resolution  of  the  1-megacycle  clock  used  in  the  test 
setup. 

-  The  normalized  stability  error  with  respect  to  time  and  shutdown  was  approxi¬ 
mately  1  x  10"4. 

-  The  quick  start  tests  indicated  that  the  RPA  essentially  reached  its  steady-state 
indication  in  about  5  minutes. 

-  An  analysis  of  temperature  sensitivity  indicates  that  perfect  temperature  com¬ 
pensation  is  theoretically  feasible  In  a  practical  design. 

In  general,  the  errors  measured  during  the  stability  and  quick  start  tests  are 
considered  to  be  more  representative  of  the  limitations  of  the  feasibility  model  and  its 
associated  tost  equipmont  than  of  the  capability  of  n  prototype  RPA.  It  Is  strongly 
believed  that  the  results  of  the  test  program  confirm  the  excellent  stability  and  quick 
start  characteristics  predicted  on  the  basis  of  an  error  analysis  of  the  RPA. 

6.2  RECOMMENDATIONS 

Based  upon  the  favorable  results  of  the  test  program,  it  is  recommended  that  a 
follow-on  program  be  funded  to  include  the  tasks  described  in  the  following  paragraphs. 

6.2,1  Design  and  Fabrication  of  an  Engineering  Prototype 

Inasmuch  as  having  a  variable  period  has  a  number  of  undesirable  consequences, 
including  a  substantial  computing  problem,  it  is  recommended  that  an  engineering  prototype 
RPA  be  designed  and  fabricated  to  include  a  synchronous  driving  motor  connected  to  the 
rotor  via  a  very  flexible  coupling  which  would  have  negligible  effect  upon  the  modulation 
resulting  from  an  input  acceleration.  The  stability  and  linearity  of  the  spring  constant 
should  have  negligible  effect  upon  the  accuracy  of  the  synchronous  RPA,  since  the  torque 
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applied  by  the  spring  will  be  only  sufficient  to  balance  the  damping  torque.  It  is  shown 
in  paragraph  5.2  that  substantial  changes  in  the  damping  torque  have  negligible  effect 
upon  the  accuracy  of  the  RPA.  The  additional  noise  introduced  by  the  hunting  of  the 
synchronous  motor  is  not  expected  to  affect  the  accuracy  potential  of  the  RPA  to  a 
significant  degree. 

6.2.2  Test  Program 

A  functional  and  environmental  test  program  should  be  conducted  on  the  engineering 
model  to  determine  stability  for  the  dormant  missile  application. 

6.2.3  Further  Modifications 

A  study  of  further  modifications  to  the  acceleration  sensor,  to  the  auxiliary  equip¬ 
ment,  and  to  data  processing  methods  should  be  made.  These  modifications  would  result 
in  improved  performance  and  accuracy,  shorter  reaction  time  or  equipment  simplifica¬ 
tions.  This  study  (intended  to  be  performed  in  parallel  with  the  design  and  test  programs) 
would  result  in  a  report  giving  analytical  data  and  recommendations  for  further  modifi¬ 
cations. 


